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Abstract. We present a reanalysis of the light curves of HD 23642, a detached eclipsing binary star in the Pleiades open cluster, 
with emphasis on a detailed error analysis. We compare the masses and radii of the two stars to predictions of stellar evolutionary 
models and find that the metal and helium abundances of the Pleiades are approximately solar. We present a new method for 
finding distances to eclipsing binaries, of spectral types A to M, using the empirical calibrations of effective temperature versus 
surface brightness given by Kervella et al. 1 2004 1. We use the calibration for A*-filter surface brightness to determine a distance 
of 139.1 + 3.5 pc to HD 23642 and the Pleiades. This distance is in excellent agreement with distances found from the use 
of theoretical and empirical bolometric corrections. We show that the determination of distance, both from the use of surface 
brightness relations and from the use of bolometric corrections, is more accurate and precise at infrared wavelengths than at 
optical wavelengths. The distance to HD 23642 is consistent with that derived from photometric methods and Hubble Space 
Telecscope parallaxes, but is inconsistent with the distance measured using Hipparcos parallaxes of HD 23642 and of other 
Pleiades stars. 

Key words, stars: distances - stars: fundamental parameters - stars: binaries: eclipsing - open clusters and associations: general 
- stars: binaries: spectroscopic - stars: chemically peculiar 



1. Introduction 

The Pleiades is a nearby, young open star cluster which is of 
fundamental importance to our understanding of stellar evo- 
lution and the cosmic distance scale. It has been exhaustively 
studied by many researchers and its distance and chemical 
composition were, until recent observations, considered to be 
well established. The distance derived from data obtained by 
the Hipparcos satellite, however, is in disagreement with tradi- 
tional values, leading to claims that stellar evolutionary theory 
is much less reliable than previously thought. 

The 'long' distance scale of 132 + 3 pc was established by 
main sequence fitting analyses (e.g., Johnson 1957 Meynet, 
Mermilliod & Maeder 1993). Recent parallax observations 
from terrestrial telescopes (Gatewood, de Jonge & Han 2000 1, 
and from the Hubble Space Telescope (Benedict et al. [2004) 
are in good agreement with this distance. 

The astrometric binary HD 23850 (Atlas) was recently 
studied by Pan, Shao & Kulkarni (2004) using the Palomar 
Testbed Interferometer (Colavita et al. 1999 1. These authors did 
not have a spectroscopic orbit for HD 23850 available to them, 
but were able to show that the distance to Atlas was greater than 
127 pc, and probably between 133 and 137 pc. Zwahlen et al. 



new interferometric measurements which, combined with the 
observations of Pan et al. |2004), give an entirely geometrical 
distance of 132 + 4 pc to HD 23850. 

A 'short' distance scale of 120 + 3 pc (van Leeuwen 2004) 
has been found using trigonometrical parallaxes observed by 
the Hipparcos space satellite (Perryman et al. 1997). This is 
in conflict with the traditional 'long' distance scale for the 
Pleiades. In an attempt to explain this, vanLeeuwen (1999) 
placed the main sequences of other nearby open clusters in the 
HR diagram using Hipparcos parallaxes, and found that five of 
the eight clusters have main sequences as faint as the Pleiades. 

Castellani et al. (2002) have shown that current theoreti- 
cal stellar evolutionary models can fit the Pleiades main se- 
quence if a low metal abundance of Z — 0.012 is adopted. 
However, Stello & Nissen (2001 ) used a metallicity-insensitive 
photometric technique to demonstrate that, if the Hipparcos 
parallaxes were correct, the main sequence Pleiades stars were 
implausibly fainter than their counterparts in the field. Also, 
Boesgaard & Friel ( 1990) have measured the iron abundance 
of the Pleiades to be approximately solar ([Fe/H] = -0.034 + 
0.024) from high-resolution spectra of twelve F dwarfs in the 
cluster (further references can be found in Stauffer et al. 2003 ). 



2004 ) have subsequently published a spectroscopic orbit and 
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Narayanan & Gould (1999) have presented evidence that 
the Hipparcos parallaxes are correlated on angular scales of 
two to three degrees. They used a variant of the moving 
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cluster method, based on radial velocities, to find a distance 
of 130 + llpc, in agreement with both the 'long' distance 
scale and the 'short' Hipparcos distance (van Leeuwen l2004> . 
Makarov (2002) has reanalysed the Hipparcos data, allowing 
for this suggested correlation, and found the Pleiades distance 
to be 129 + 3pc. Until this result is confirmed, however, the 
'long' and 'short' distance scales cannot yet be considered to 
be reconciled. 

Munari et al. (2004 hereafter M04) studied the detached 
eclipsing binary HD 23642 and found a distance of 132 + 2 pc, 
in good agreement with the 'long' distance scale. The method 
used by M04 is commonly used to find the distances to eclips- 
ing binaries but depends on theoretical calculations to provide 
bolometric corrections. We have reanalysed the data of M04 
(which these authors have graciously provided on the internet) 
to investigate alternative, empirical, methods of finding the dis- 
tance to HD 23642 and similar eclipsing binaries by the use of 
surface brightness relations. 

1.1. The eclipsing binary HD 23642 in the Pleiades 

Detached eclipsing binaries (dEBs) with double-lined spectra 
are one of the best sources of fundamental astrophysical data 
(Andersen ll991> because their absolute masses and radii can be 
measured to accuracies better than 1%. Such data can be used 
to provide a strict test of different stellar evolutionary models, 
and the distances to dEBs can be determined empirically to an 
accuracy of about 5%. dEBs in open clusters are particularly 
useful because the age and chemical composition of the cluster 
can be combined with accurate values of the masses and radii 
of the dEB to provide an even more exacting test of theoretical 
models (Southworth, Maxted & Smallev l2004all2004b> . 

HD 23642 (Table ^ was discovered to be a double-lined 
spectroscopic binary by Pearce fT9~57i and Abt (1958 1, and 
both components have been found to display slight spectral pe- 
culiarities (Abt & Levato [T978> . Torres ( 2003 ) discovered shal- 
low secondary eclipses in the Hipparcos photometric data of 
HD 23642 and also presented an accurate spectroscopic orbit. 
M04 derived precise absolute masses and radii of both com- 
ponents from high-resolution spectra and complete BV light 
curves. M04 found a distance of 131.9 + 2.1 pc, in disagree- 
ment with the Hipparcos parallax distance of 111 + 12pc for 
HD 23642. 



Table 1. Identifications and astrophysical data for HD 23642. 
References: (1) Perryman et al. dl997> : (2) Abt & Levato 
(I1978i : (3) 2 Micron All Sky Survey; (4) M04; quantities in 
parentheses are uncertainties in the final digit of the quanities. 





HD 23642 


References 


Hipparcos number 


HIP 17704 


1 


Hipparcos distance (pc) 


111 ± 12 


1 


Spectral type 


AOVp (Si) + Am 


2 


Bj 


6.923 ± 0.011 


1 


V T 


6.839 ±0.011 


1 


^2MASS 


6.635 ± 0.023 


3 


#2MASS 


6.641 ± 0.026 


3 


^2MASS 


6.607 ± 0.024 


3 


Orbital period (days) 


2.46113400(34) 


4 


Reference time (HJD) 


2452903.5981(13) 


4 




0.0 0.2 0.4 0.6 0.8 1.0 

Phase 



Fig. 1. Spectroscopic orbit for HD 23642 



Table 2. Spectroscopic orbital parameters for HD 23642. 



Primary Secondary 
Semiamplitude K (kms~') 99. 10 ±0.58 140.20 ±0.57 
Systemic velocity ( km s~ 1 ) 6.07 ± 0.39 

Mass ratio q 0.7068 ± 0.0050 

asin/(R G ) 11.636 ±0.040 

M sin 3 i ( M G ) 2.047 ±0.021 1.447 ±0.017 



2. Spectroscopic analysis 

M04 observed HD 23642 five times with the Elodie echelle 
spectrograph on the 1 .93 m telescope of the Observatoire de 
Haute-Provence. The radial velocities derived were combined 
with the spectroscopic observations of Pearce ( 1957 1 and Abt 
( 1958 1, using lower weights for the older data, to calculate a 
circular spectroscopic orbit. 

The low weight - and low precison - of the data of Pearce 
( 119571 and Abt d!958> mean that they contribute little to the 
accuracy of the spectroscopic orbit. For comparison with the 
results of M04 we have chosen to derive the orbit using only 
the five echelle velocities for each star. The orbit was com- 



puted using sbop 1 , with the orbital ephemeris from M04, ec- 
centricity fixed at zero, and equal systemic velocities for both 
stars. The root-mean-squares of the residuals of the resulting 
spectroscopic orbit are 0.4 and 1.2kms _1 for the primary and 
secondary stars, respectively. The spectroscopic orbit is plotted 
in Fig. ^and its parameters are given in Table |2] The orbital 
parameters are in acceptable agreement with those of M04 and 
Torres ( 120031 . 



1 Spectroscopic Binary Orbit Program written by Dr. P. B. Etzel 
(http : //mintaka . sdsu.edu/faculty/etzel/). 
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Fig. 2. Comparison between a spectrum of HD 23642 and the best-fitting synthetic spectrum used to determine the atmospheric 
parameters of the stars. The spectrum of HD 23642 plotted here is a recombination of the individual spectra of the two stars, which 
were obtained by spectral disentangling (Simon & Sturm 1994). The source ion of some lines of interest have been indicated, 
with the rest wavelength of the line (A) and which star is producing it (A for the primary or B for the secondary). The effective 
radial velocities of the primary and secondary stars in this diagram are -58 and +160kms~ 1 , respectively. 



2.1. Determination of effective temperatures 

Atmospheric parameters were derived for the components of 
HD 23642 by comparing the observed spectra (from M04) with 
synthetic spectra calculated using uclsyn (Smith 1992 Smalley 
et al. 2001 1, Kurucz (1993 1 atlas9 model atmospheres with- 
out approximate convective overshooting (Castelli, Gratton & 
Kurucz 1997) and absorption lines from the Kurucz & Bell 
( 1995 ) linelist. The spectra were rotationally broadened as nec- 
essary and instrumental broadening was applied with a FWHM 
of 0.11 A to match the resolution of the observations. Surface 
gravities of 4.25 were assumed for both stars. 

For the primary star, spectroscopic fitting gives an effec- 
tive temperature T e s,K — 9750 + 250 K with a microturbu- 
lent velocity of £t,a = 2kms~ 1 and a rotational velocity of 
Va sin i = 37 ±2 km s _I . For the secondary star we find r e tf,B 
7600 + 400 K, ^t,b = 4kms M and V B sin; = 32 + Skms -1 . 
The values for microturbulence velocity are consistent with 
those typically found for stars of these effective temperatures 
(Smalley 2004). The quoted uncertainties are limits of high 
confidence and are somewhat larger than the formal fitting er- 
rors. A monochromatic light ratio of 0.25 + 0.05 was obtained 
at 4480 A. The observations and the best-fitting synthetic spec- 
trum are represented in Fig.0 

Atmospheric parameters can be also estimated from pho- 
tometric indices. We have obtained uvby/3 photometry from 
Hauck & Mermilliod ( 1988) and dereddened it using E^-y = 
0.008, calculated from E B _ V = 0.012 (M04) and E b _ y * 
0.73E B _ V (Crawford ll975> . Using the semi-empirical grid cal- 
ibrations of Moon & Dworetsky ( 1985 ) and the tefflogg pro- 
gram (Moon [T985l . we obtained T eS = 9200 and log g = 4.30, 
for the combined light of the system. To evaluate the effects 
of the secondary we have subtracted the photometry of the 
classical Am star 63Tauri using a V-filter magnitude differ- 



ence of 1.44. This gave the parameters T e gj^ - 9870 K and 
log gA - 4.37 for the primary component, in good agreement 
with our observationally determined parameters for this star. 

A near-fundamental determination of effective temperature 
can be obtained using the Infrared Flux Method (Blackwell & 
Shallis I1977> . Ultraviolet fluxes were obtained from the IUE 
archive, optical fluxes from Kharitonov et al. (1988) and in- 
frared fluxes from the 2MASS catalogue. From this the total 
integrated flux at the Earth was found to be (5.44 + 0.44) x 10~ 8 
ergs -1 cirT 2 , for a reddening of Eb-v = 0.012 (taken from 
M04). The IRFM then yielded 77 eff = 8900 + 350 K, which 
is rather low compared to the above values but is affected by 
the flux contribution of the cooler secondary star, which is pro- 
portionately brighter in the infrared. Allowing for the presence 
of the secondary star using the method of Smalley ( 1993} we 
find that the primary would have r e ff,A = 9250 + 400 K for a 
secondary star with r e ff,B = 7500 + 500 K, which is consistent 
with the values determined above. 

Fundamental effective temperatures can be obtained for bi- 
nary systems using total integrated fluxes and angular diam- 
eters obtained from system parameters, and known distances 
(Smalley & Dworetsky [1995J Smalley et al. l2002> . In the case 
of HD 23642 the properties of the system have been found us- 
ing a model-dependent method, so application of this proce- 
dure would lead to a circular argument. However, the method 
does allow for a consistency check on the two effective temper- 
atures and, importantly, their error estimates. Using the param- 
eters obtained in the present work, we find r e ft,A = 9620 + 
280 K and 77 eff3 = 7510 + 430 K for the primary and sec- 
ondary, respectively. Similar results are obtained for the param- 
eters given by M04. However, use of the Hipparcos parallax of 
HD 23642 (which gives a distance of 1 1 1 + 12 pc) would give 
r effjA = 8640 + 540 K and 77 eff , B = 6690 + 570 K, which are 
clearly inconsistent with the values obtained above. The 'short' 



4 



J. Southworth et al.: HD 23642 and the distance to the Pleiades 



Pleiades distance ( 1 20 + 3 pc) would give r effjA = 9000 + 310K 
and r e ff,B = 6970 ± 450 K, which is closer but still somewhat 
discrepant. 

Using several techniques we have found the effective tem- 
peratures of the two components stars of HD 23642 to be 
Tdf.A = 9750 + 250 K for the primary and r eff , B = 7600 + 400 K 
for the secondary. Our error estimates are higher than those re- 
ported in M04, primarily because we have assessed the influ- 
ence of external uncertainties, in addition to the internal preci- 
sion of fits to spectra. 



3. Photometric analysis 

The B and V light curves contain 432 and 492 individual 
measurements, respectively, obtained with a 28 cm Schmidt- 
Cassegrain telescope and photometer by M04. The two light 
curves were solved separately using ebop (Nelson & Davis 
119721 Popper & Etzel 1981 1. This is a simple and efficient light 
curve fitting code in which the discs of the stars are modelled 
using biaxial ellipsoids. Linear limb darkening coefficient val- 
ues of 0.496 and 0.596 (B) and 0.421 and 0.548 ( V), for the pri- 
mary and secondary stars respectively, were adopted from van 
Hamme ( 1993 } as the light curves are not of sufficient quality 
to include them as free parameters. Gravity darkening expo- 
nents B\ were fixed at 1.0 (Claret 1998) and the mass ratio was 
fixed at the spectroscopic value. The ephemeris given in M04 
was used and the orbit was assumed to be circular. 

Although the M04 light curves are of reasonable quality, 
deriving accurate parameters from them is problematic due to 
the shallow eclipses. This is exacerbated by some scattered data 
in the B light curve, which makes it less reliable than the V 
light curve. As contaminating 'third' light, Lj, is poorly con- 
strained by the observations, we have made separate solutions 
for L3 = and 0.05 (in units of the total light of the eclipsing 
stars) and included differences in the parameter values derived 
in the uncertainties quote below. As there are no features in 
the spectra of HD 23642 known to come from a third star, it is 
unlikely that third light is greater than 5%. 

Initial solutions provided an inadequate fit to the light vari- 
ation outside eclipse so the reflection effect for the secondary 
star was separately adjusted towards best fit rather than being 
calculated from the system geometry. We also solved the light 
curves using the Wilson-Devinney code (Wilson & Devinney 
[T97T1 Wilson [T993l . using the 1998 version (wd98) with a de- 
tailed treatment of reflection. As the differences between the 
ebop and wd98 solutions were negligible, further analysis was 
undertaken using ebop. This code has two important advan- 
tages; a detailed error analysis is not prohibitively expensive 
in terms of computer time, and the philosophy of the ebop code 
is to solve for the set of parameters most directly related to the 
light curve shape. 

As the photometry of M04 is not supplied with observa- 
tional errors, we have weighted all observations equally. We 
will judge the quality of the fit of an ebop model light curve to 
the observational data using the root mean square of the resid- 
uals of the fit, cr rms . 



3.1. Light curve solution 

The ratio of the radii, k, is poorly constrained by the light 
curves because the eclipses are very shallow. Whilst a reason- 
able photometric solution can be obtained from the light curves 
alone, an alternative is to use a spectroscopic light ratio to con- 
strain k. A light ratio of — — 0.31 + 0.03 was given by Torres 
(2003), based on a cross-correlation analysis of a 45 A wide 
spectral window centred on 5 187 A. This spectroscopic light 
ratio is noted to be preliminary so we provide separate solu- 
tions without ('solution A') and with ('solution B') its inclusion 
in the light curve fitting procedure, but preference is given to 
solutions including the spectroscopic constraint. We have used 
the light ratio of Torres (2003 ) as it is based on a larger amount 
of observational data than the light ratio found in Section ITTl 
and because the wavelength it was obtained at is closer to the 
central wavelength of the V filter. 

The light ratio of Torres (2003) was converted to a V fil- 
ter light ratio using a V filter response function and synthetic 
spectra, calculated from atlas9 model atmospheres, for the ef- 
fective temperatures and surface gravities found in our prelimi- 
nary analyses. The resulting V filter light ratio of 0.335 ± 0.035 
(where the uncertainties include a small contribution due to 
possible systematic errors from the use of atlas9 model atmo- 
spheres) has been used to constrain k using the V light curve. 
The resulting values of k were then adopted for solution of the 
B light curve. 

Table|5]gives solution A, and Table|4]gives solution B. The 
best fit for the former solution is compared to the observational 
data in Fig.|5J the light variation for the latter solution is almost 
identical so has not been plotted. Note that the B light curve 
appears to be badly fitted at the centre of the secondary eclipse. 
Investigation has revealed that this is not a problem with the 
ebop model, but is caused by scatter present in the observa- 
tional data. A better fit can be obtained by rejecting one night's 
data, around phase 0.48, which is brighter than the model light 
curve. After rejection of these data, the fit is significantly im- 
proved but the derived parameters are quite similar. We have 
therefore included all observational data in the fitting proce- 
dure, and suggest that further photometric data should be ob- 
tained. 

The uncertainties in the fitted parameters were estimated 
using Monte Carlo simulations (see Southworth et al. 2004b 
for details). In this algorithm, a best fit is found and the ob- 
servational scatter is calculated. The model is then evaluated at 
the phases of observation, random Gaussian noise (of the same 
size as the actual observational scatter) is added, and the result- 
ing synthetic light curve is fitted. This procedure is undertaken 
10 000 times for each observed light curve. The uncertainties in 
the photometric parameter values are estimated by calculating 
the standard deviation of the values found during the Monte 
Carlo simulations. The resulting uncertainties are included in 
Tableland Table 0] Some results of the Monte Carlo simu- 
lations are shown in Fig.@] Uncertainties in the theoretically- 
derived limb darkening coefficients have been incorporated by 
perturbing the values of the coefficients by +0.05, on a flat dis- 
tribution, for each Monte Carlo simulation. The uncertainties 
estimated using this Monte Carlo simulation algorithm have 
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Fig. 3. The M04 B and V light curves with our best fitted overplotted. The V light curve is shifted by +0. 1 mag for clarity. The 
residuals of the fit are offset by +0.24 mag and +0.28 mag for the B and V light curves respectively. Note that the poor fit around 
the secondary eclipse in the B light curve is due to scattered data, as suggested by the distribution of the residuals for this light 
curve. Rejection of the offending data makes the fit look better but is otherwise unjustified (see text for discussion). 

Table 3. Photometric solution A, found by fitting the light curves without using a spectroscopic light ratio to constrain the ratio 
of the radii. The individual errors are from Monte Carlo simulations. The combined parameters (bottom line) are the weighted 
means of the two different values for each parameter. The quality of the fit is given by cr rms (see text for definition). 



Light 


Surface bright- 


Light 


Primary radius, 


Ratio of the 


Secondary radius, 


Orbital 


0" rms 


curve 


ness ratio, J 


ratio 


r\ (a) 


radii, k 


rs (a) 


inclination, i (°) 


(mmag) 


B 


0.316 + 0.018 


0.264 + 0.083 


0.1536 ±0.0069 


0.929 ±0.1 16 


0.1427 ±0.0108 


76.91 ±0.41 


4.639 


V 


0.493 ±0.018 


0.366 + 0.094 


0.1487 ±0.0057 


0.882 ± 0.093 


0.13 12 ±0.0083 


77.80 ±0.24 


3.254 


Combined 






0.1507 ±0.0044 


0.900 ±0.073 


0.1355 ±0.0066 


77.57 ±0.21 





been found by the authors to be extremely reliable (Southworth 
et al. l2004cll2"nolbl . 

The fractional stellar radii given by solution A are r& = 
0.1507 ± 0.0044 and r B = 0.1355 ± 0.0066, whereas inclusion 
of the spectroscopic light ratio (solution B) gives r& = 0.1538 ± 
0.0024 and r B = 0.1 300 ±0.0037. The fractional radii found by 
M04werer A = 0.1514 + 0.0025 and m = 0.1254 + 0.0022. Our 
spectroscopic-constraint result agrees well with the M04 result, 
but we are unable to reproduce the small uncertainties claimed 
by M04. Fig.[5]compares the values of k found for solution A, 
solution B and by M04, to the residuals of the fit for a range 
of different values of k for each light curve. Solution A has the 
largest uncertainty but is close to the minima of the residuals 
in the B and V light curves. Solution B, for which k was found 
using a spectroscopic light ratio, has a smaller uncertainty in k 
but has a slight dependence on theoretical model atmospheres. 
The results of M04 were obtained from the same data as our 
solution A, but the two values of k are quited different. The low 



uncertainties quoted by M04 mean that their value of k (calcu- 
lated by us from the stellar radii given by M04) is inconsistent 
with the minima in the residuals curves. M04 adopted formal 
errors from their photometric analysis, which are known to be 
generally somewhat optimistic (Popper 1984 1. 

4. Absolute dimensions and comparison with 
stellar models 

The absolute dimensions of the two stars have been calculated 
from the spectroscopic results and photometric solution B, and 
are given in Table |5] We note that the radii of the two stars 
calculated using solution A are = 1.796 ± 0.053 R Q and 
Rb = 1.615 + 0.079 R Q . The masses and radii of the two com- 
ponents of HD 23642 can be compared to stellar evolutionary 
models to estimate the metal abundance of the system. This is 
important because unusual chemical compositions have been 
suggested as possible reasons why the 'long' distance scale of 
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Table 4. Photometric solution B, found using ebop and the light ratio of Torres (2003 1. The individual uncertainties are from 
Monte Carlo simulations with a fixed ratio of the radii, k. Symbols have the same meaning as in Table[3] The adopted parameter 
values (bottom line) are the weighted means of the different parameter values for each light curve. The adopted uncertainties 
depend mainly on the possible values of k. 



Light 


Light 






rA 




i 


0"rms 


curve 


ratio 


k 


J 


(a) 


(a) 


o 


(mmag) 


V 


0.300 


0.8075 


0.483 ±0.014 


0.1531 ±0.0009 


0.1237 ±0.0007 


78.04 ±0.08 


3.254 


V 


0.335 


0.8477 


0.489 ±0.015 


0.1508 ±0.0009 


0.1278 ±0.0008 


77.90 ±0.08 


3.254 


V 


0.370 


0.8858 


0.493 ±0.015 


0.1485 ±0.0009 


0.1315 ±0.0008 


77.79 ±0.08 


3.254 


B 


0.1 90 ±0.008 


0.8075 


0.304 ±0.014 


0.1599 ±0.0019 


0.1291 ±0.0015 


77.47 ±0.19 


4.654 


B 


0.214 ±0.008 


0.8477 


0.309 ±0.015 


0.1580 ±0.0019 


0.1339 ±0.0016 


77.26 ±0.19 


4.647 


B 


0.237 ±0.009 


0.8858 


0.313 ±0.015 


0.1560 ±0.0019 


0.1382 ±0.0016 


77.08 ±0.19 


4.642 


Adopted 




0.848 ±0.039 




0.1538 ±0.0024 


0.1300 ±0.0037 


77.78 ±0.17 
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r(B) / r(A) 




0.130 0.145 0.160 
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0.8 1.0 
r(B) / r(A) 
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Fig. 4. Results of the Monte Carlo analysis for the B (left pan- 
els) and V (right panels) light curves. The top two panels show 
ta plotted against r%. The middle two panels show k versus the 
light ratios. The lower two panels show the reduced^ 2 versus 
k. These values were calculated using the residuals of the best 
fits to the observed light curves. Note the greater scatter of the 
Monte Carlo solutions for the B light curve, which is due to 
the larger observational errors. Only 2000 of the 10000 points 
have been plotted in each panel. 



4.67 
" 4.66 

E 
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3.28 
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-e- 



sol. A 



M04 



0.70 0.80 0.90 1.00 1.10 

r(B) / r(A) 

Fig. 5. Comparison between the quality of the fit for different 
values of k and the solutions derived in this work and by M04. 
<Xmis has been plotted against k for the B light curve (middle 
panel) and V light curve (lower panel) using the same scales. 
The optimum values of k found in the individual solutions of 
the two light curves are indicated by filled squares. The upper 
panel shows the values and uncertainties of k found in solution 
A, solution B, and by M04. Note that the quality of the pho- 
tometry is excellent, particularly for the V light curve where 
the photometric error of an individual datapoint is 3.25 milli- 
magnitudes. 



the Pleiades is in disagreement with the Hipparcos parallax dis- 
tances. 



J. Southworth et al.: HD 23642 and the distance to the Pleiades 



7 



Table 5. Absolute dimensions of the components of HD 23642, 
calculated using photometric solution B. Symbols have their 
usual meanings and the equatorial rotational velocities are de- 
noted by V eq . The absolute bolometric magnitudes have been 
calculated using a solar luminosity of 3.826xl0 26 W and abso- 
lute bolometric magnitude of 4.75. 



Primary star Secondary star 



Mass ( Mq) 


2. 193 ±0.022 


1.550 


±0.018 


Radius (Rq) 


1.831 ±0.029 


1.548 


± 0.044 


log 8 (cm s -1 ) 


4.254 ±0.014 


4.249 


±0.025 


T eS (K) 


9750 ± 250 


7600 


±400 


log(L/L s ) 


1.437 ±0.047 


0.858 


±0.095 




1.16 ±0.12 


2.60 


±0.24 


VeqdcmS-') 


38 ±2 


33 


±3 


Vsynch (kms" 1 ) 


37.6 ±0.6 


31.8 


±0.9 



In Fig. [6] the masses and radii of the components of 
HD 23642, found from solution B, have been compared to 
predictions of the Granada stellar evolutionary models (Claret 
fT995llT997l Claret & Gimenez lT995l An age of 12 5 Myr has 
been adopted (Stauffer, Schultz & Kirkpatrick 1998 1; a change 
in this age by -25 or +75 Myr does not affect the conclu- 
sions below. Fig. [6] shows predictions for metal abundances 
of Z = 0.01, 0.02 and 0.03. For each metal abundance we 
have plotted predictions for normal helium abundance (dashed 
lines) and for significantly enhanced helium abundances (dot- 
ted lines). Fig. 0compares the masses and radii of HD 23642 
to the predictions of the Cambridge stellar evolutionary models 
(Pols et al. [T998l for metal abundances of Z = 0.01, 0.02 and 
0.03. The solar chemical composition isochrone is also shown 
for an age of 175 Myr. 

The masses and radii of the components of HD 23642 sug- 
gest that the metal abundance of the dEB is the slightly greater 
than solar (Z « 0.02). Predictions for enhanced helium abun- 
dances are ruled out as they predict a mass-radius relation 
for HD 23642 much steeper than observed. The approximately 
solar Pleiades iron abundance found by Boesgaard & Friel 
( 1990 1, from high-resolution spectroscopy of F dwarfs, is con- 
firmed. The 'short' and 'long' Pleiades distances therefore can- 
not be reconciled by adopting an unusual chemical composition 
for the cluster. 

5. The distance to HD 23642 and the Pleiades 

We will now derive the distance to HD 23642 using three meth- 
ods, one of which is introduced here. We investigate the normal 
technique using bolometric corrections to find absolute visual 
magnitudes, the use of surface brightness calibrations in terms 
of observed colour indices, and a new method based on infrared 
surface brightness calibrations in terms of effective tempera- 
ture. All the techniques require reliable apparent magnitudes, 
which for HD 23642 are available from the Tycho experiment 
on the Hipparcos satellite and from 2MASS. The apparent B 
and V magnitudes observed by Tycho (Perryman et al. 119971 
have been converted to the Johnson photometric system using 
the calibration of Bessell (2000 1. The JHK photometry from 
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Fig. 6. Comparison between the observed properties of 
HD 23642 and the Granada stellar evolutionary models for 
metal abundances of Z = 0.01 (circles), Z = 0.02 (triangles) 
and Z = 0.03 (squares). Predictions for normal helium abun- 
dances are plotted with dashed lines and helium-rich model 
predictions are plotted using dotted lines. An age of 125 Myr 
was assumed. 



2MASS has been converted to the SAAO near-infrared photo- 
metric system using the calibration of Carpenter ( 2001 1. 
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Fig. 7. Comparison between the observed properties of 
HD 23642 and the Cambridge stellar evolutionary models for 
metal abundances of Z = 0.01 (circles), Z = 0.02 (triangles) 
and Z = 0.03 (squares). An age of 125 Myr was assumed 
(dashed lines), but predictions for a solar chemical composi- 
tion and an age of 175 Myr are also shown (dotted line). 
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5.1. Distance from the use of bolometric corrections 

The traditional method of determining the distance of an eclips- 
ing binary is to calculate the luminosity of each star from its ra- 
dius and effective temperature. The resulting values of absolute 
bolometric magnitude, Mboi, are then converted to absolute vi- 
sual magnitudes, My, using bolometric corrections (BCs). The 
combined My of the two stars is then compared to the appar- 
ent visual magnitude, my, to find the distance modulus. This 
method is quite capable of yielding accurate results provided 
that care is taken to use effective temperature measurements 
consistent with the fundamental definition of effective temper- 
ature. Consistent zeropoints must also be used for the calcu- 
lation of absolute bolometric magnitude and the BCs (Bessell, 
Castelli & Plez [T998l . 

To find the distance to HD 23642, we have adopted the 
astrophysical parameters of the system given in Table [5] An 
interstellar reddening of Eb-v = 0.012 + 0.004 mag has 
been adopted from M04. We have calculated the distance 
to HD 23642 using the empirical BCs given by Code et al. 
dl976> . with a calibration uncertainty of 0.05 mag (Clausen 
2004) added in quadrature. We have also derived the distance 
using the theoretically-calculated BCs of Bessell et al. (1998 1 
and of Girardi et al. (2002 1. For the primary component of 
HD 23642, the uncertainty resulting from its effective tempera- 
ture measurement is reduced due to the form of the BC function 
around 10000K. The main contribution to the overall uncer- 
tainty comes from the uncertainties in the effective tempera- 
tures of the stars. 

The distances found by using BCs are given in Table |6] 
and show that a value around 139pc is obtained consistently. 
The distance found using the empirically derived BCs of Code 
et al. ( 1976 1 is very similar to the distances found by using 
the theoretically-derived BCs, suggesting that systematic errors 
due to the use of model atmospheres are small. Distances de- 
rived using BCs for the JHK filters are more precise because 
the uncertainties in effective temperature and interstellar red- 
dening are smaller. If the results of photometric solution A are 
used to find the distance to HD 23642 using this method, dis- 
tances of around 139pc are also found, but with slightly larger 
uncertainties. 

The BCs given by Girardi et al. (2002 1 are available for sev- 
eral different metal abundances, [M/H]. We have investigated 
the effect of non-solar metal abundances on the distance de- 
rived using BCs for [M/H] = -0.5 and +0.5 (Table|6j. Whilst 
a significant change in distance is found for distances derived 
using BCs for the V filter, the effect is much smaller for the 
K filter, underlining the usefulness of infrared photometry in 
determining distances to dEBs. 

M04 found the distance to HD 23642, using the Bessell et 
al. ([1998) BCs, to be 131.9+2.1 pc. However, we have been un- 
able to reproduce their result using the properties of HD 23642 
found by these authors. The absolute bolometric magnitudes 
given by M04 appear to have been taken from output files 
produced by the Wilson-Devinney light curve modelling code, 
which uses a solar luminosity of L = 3.906 x 10 26 W and 
absolute bolometric magnitude of Mboi© = 4.77. The BCs 
used by M04, however, were calculated using different values: 



Table 6. The distances derived for HD 23642 by using several 
different sources of bolometric corrections. 



Bolometric corrections 


Filter 


[M/H] 


Distance (pc) 


Code et al. <1976> 


V 




138.1 


+ 6.2 


Bessell et al. fl998l 


V 




uy. j 




Bessell et al. Jl998l 


K 




1 JO.U 


X J.J 


Girardi et al. t2002> 


B 


0.0 


140.2 


+ 6.1 


Girardi et al. i2002l 


V 


0.0 


139.8 


±5.3 


Girardi et al. i2002l 


J 


0.0 


138.7 


±3.8 


Girardi et al. <2002l 


H 


0.0 


138.0 


±3.3 


Girardi et al. ( 2002 1 


K 


0.0 


138.8 


+ 3.3 


Girardi et al. ( 2002 1 


V 


-0.5 


137.9 


±5.3 


Girardi et al. i2002l 


V 


+0.5 


142.4 


±5.3 


Girardi et al. i2002l 


K 


-0.5 


137.9 


±3.3 


Girardi et al. {20021 


K 


+0.5 


139.8 


±3.3 



L = 3.855 x 10 26 W and M bol0 = 4.74 (Bessell et al. 1998 1. 
This inconsistency appears to be sufficient to explain the dis- 
crepancy between their result and our results. If we adopt the 
masses, radii and effective temperatures of M04, we find a dis- 
tance of 135.5 ± 2.3 pc. This is about 3pc smaller than most 
of our results, mainly due to the smaller stellar radii found by 
M04 for the components of HD 23642. 

The discussion above shows that it is possible to calculate 
consistent distances to eclipsing binaries, using different sets of 
BCs. One weakness of this method is that it is difficult to eval- 
uate the systematic error introduced by the uncertainty in the 
zeropoint of the effective temperature scale. The use of BCs 
derived from theoretical models also introduces a systematic 
error which is likely to be small in this case but, in general, 
is difficult to quantify. This systematic error is due to deficien- 
cies in the model, e.g., the approximate treatment of convection 
and the lack of complete spectral line lists. It may also be the 
case that the star has properties which are not accounted for 
by the model, e.g., spectral peculiarity due to magnetic fields 
or slow rotation. For these reasons it is desirable to develop an 
empirical distance determination method which is less sensi- 
tive to systematic errors in effective temperature and in which 
the sources of uncertainty are more explicit. 

5.2. Distance from the use of relations between 
surface brightness and colour indices 

The concept of surface brightness was first used in the anal- 
ysis of dEBs about one hundred years ago (Kruszewski & 
Semeniuk ll999> . Effective temperatures of the components of 
several dEBs were found from the stellar radii and distances de- 
termined using parallax measurements, and relations between 
surface brightness and spectral type were investigated. 

Barnes & Evans i 19761 Barnes, Evans & Parsons (1976 1 
provided a calibration between visual surface brightness, Sy, 
and the V - R colour index. The Barnes-Evans relation was 
used by Lacy J1977I to find the distances to nine dEBs. 
Popper ( 1980 1 presented calibrations between S y and the B-V 
and b — y photometric indices. Surface brightness calibrations 
for broad-band filter indices have recently been presented by 
Fouque & Gieren ( 1997 1, Di Benedetto ( 1998 1, Beuermann et 
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al. i fTgggt and Ker vellae t al. (2004 hereafter KTDS04). Salaris 
& Groenewegen (2002) have calibrated the relation between 
S v and the Stromgren photometric index, C\ , for determination 
of the distances to B stars in the Magellanic Clouds. 

In this method the angular diameter of the star is calculated 
and compared to its linear diameter, determined from photo- 
metric analysis, to find the distance. The advantage of this pro- 
cedure is that the calculations can be entirely empirical, de- 
pending on the way in which the surface brightness calibra- 
tion is constructed. One disadvantage is that individual appar- 
ent magnitudes and colour indices of the components of the 
dEB, calculated from the light ratios found in light curve anal- 
yses, must be used. The uncertainties in these quantities can 
cause the derived distance to have a low precision. 

We have determined a distance to HD 23642 of 138 + 19pc 
using the S y versus B - V calibration given by Di Benedetto 
dl998> . The main uncertainty in this result comes from the un- 
certainties in the light ratios in the B and V filters. The most 
recent calibration, given by KTDS04, has not been used be- 
cause it does not allow for the nonlinear dependence of S y on 
B-V. The B - V colour index is also not a good indicator of sur- 
face brightness because the B filter is known to be sensitive to 
metallicity through the effect of line blanketing. The V - K and 
B - L indices are good surface brightness indicators (KTDS04; 
Di Benedetto 1998 1 because they are more sensitive to surface 
brightness and because the intrinsic scatter in the calibrations 
falls below 1%. 

5.3. Distance from the use of relations between 
surface brightness and effective temperature 

Empirical relations between effective temperature and the 
wavelength-dependent surface brightness of a star, S m where 
mi is an apparent magnitude in filter A, have been derived by 
KTDS04 from interferometric observations. This allows the 
surface brightnesses of the components of a dEB to be found 
without using the light ratios of the system found during the 
light curve analysis. 

The zeroth-magnitude angular diameter is a measure of sur- 
face brightness and is defined to be the angular diameter a star 
would have if nix - 0: 



A ('«i=0) 



10^ 



(1) 



where <p is the actual angular diameter of a star (van Belle 
1999 1. From consideration of the definitions of surface bright- 
ness and angular diameter, it can be shown that the distance, d, 
to a dEB is given by 



d= 10' 



0.2mi 



2R> 



Am x =Q) 
J A 



2R P 



(2) 



=0) 



are in arc- 



where the stellar radii and Rb are in AU, i 
seconds and distance is in parsecs. 

Calibrations for ^- mi=0 > are given in terms of effective tem- 
perature by KTDS (where they are denoted using ZMLD^). 
These calibrations are for the broad-band UBVRIJHKL filters 



Table 7. The results and individual error budgets for distance 
estimates using the effective temperatures and overall apparent 
magnitudes of the HD 23642 system in the BVJHK filters. All 
distances are given in parsecs and the total uncertainties are the 
sums of the individual uncertainties added in quadrature. 



Uncertainty source 


b 


V 




u 

tl 


K 


spectroscopic K A 


0.3 


0.3 


0.3 


0.3 


0.3 


Spectroscopic 




U, J 


U, J 


U, J 


U, J 


Orbital inclination, i 


0.1 


0.1 


0.1 


0.1 


0.1 


Fractional radius, r& 


1.8 


1.7 


1.5 


1.4 


1.4 


Fractional radius, 


0.8 


1.0 


1.3 


1.4 


1.5 


Primary 


4.8 


3.3 


1.7 


0.7 


0.7 


Secondary T c (f 


3.8 


3.5 


2.1 


1.5 


1.4 


Reddening E B -v 


1.1 


0.8 


0.3 


0.2 


0.1 


Apparent magnitude 


1.0 


1.0 


1.9 


1.9 


1.9 


"Cosmic" scatter 


9.0 


8.3 


1.7 


1.9 


1.4 


Total uncertainty 


11.1 


9.9 


4.2 


3.8 


3.5 


Distance 


142.8 


141.4 


139.6 


138.4 


139.1 



and are valid for effective temperatures between 10000K and 
3600 K for main sequence stars. The JHKL filter calibrations 
have the least scatter and are also the least affected by interstel- 
lar reddening. For the K filter the scatter around the calibration 
is undetectable at a level of 1% so we conservatively adopt a 
scatter of 1%. We have applied the calibrations for B and V 
(derived from Tycho data) and for J, H and K (derived from 
2MASS data). As there is no "standard" infrared photometric 
system, the calibrations of KTDS04 use data from several dif- 
ferent JHKL systems, so the systematic uncertainty of not hav- 
ing a standard system is already included in the quoted scatter 
in the calibrations. For the A stars only, the scatter around the 
B and V calibrations is much smaller than the overall scatter 
quoted, so for HD 23642 the B and V distance uncertainties are 
overestimated by a factor of about two. 

The distances found using equation|2]and the KTDS04 cal- 
ibrations are given in TableQ We have calculated the distances 
from the results of the light curve solution and the spectro- 
scopic velocity semiamplitudes (K& and K%), in order to care- 
fully assess how important the uncertainties in each quantity 
are to the final distance uncertainty. The intrinsic scatter in the 
calibrations is marked as "cosmic" scatter in Tableland is the 
main contributor to the uncertainty in the B and V filter dis- 
tances. Note that the uncertainties in the JHK calibration dis- 
tances are much smaller than in the BV calibration distances, 
because the calibrations have much smaller scatter and redden- 
ing is less important. We will adopt the K filter calibration 
distance of 139.1 + 3.5 pc as our final distance to HD 23642 
and so to the Pleiades. Note that we cannot treat any of the 
distance estimates investigated above as being independent of 
each other as they are all calculated using the same values for 
reddening, stellar radii and effective temperatures. If we adopt 
the results of photometric solution A, the K filter distance we 
find is 139.7 pc with an error of 4.7 pc, which is entirely con- 
sistent with our adopted distance of 139.1 + 3.5 pc. 

One shortcoming of finding distances using equation [2] is 
that the effective temperature scales used in analysis of the dEB 
and for the calibration must be the same to avoid systematic 
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Table 8. Summary of the different distances found for the 
Pleiades or for Pleiades members, both from the literature (up- 
per part of the table) and for HD 23642 in this work (lower part 
of the table). References are also given in the text and abbrevi- 
ations and symbols have their usual meanings. 



Source of distance measurement 


Distance (pc) 


Hipparcos parallaxes (van Leeuwen 2004 1 


120 d 


: 3 


Hipparcos parallax ot HD 23642 


111 d 


: 12 


Hipparcos parallax ot HD 23850 


117 d 


: 14 


MS fitting (Stello & Nissen 2001 1 


132.4 d 


: 1.8 


Ground-based parallaxes 


130.9 i 


: 7.4 


HST parallaxes of three Pleiades stars 


134.6 d 


: 3.1 


Narayanan & Gould fHTO 


131 d 


: 11 


Makarov <2002t 


129 d 


: 3 


Eclipsing binary HD 23642 (Munari et al.|2004 1 


131.9 i 


: 2.1 


Astrometric binary HD 23850 


132 d 


: 4 


Code et al. 1 1976 1 empirical BCs (V filter) 


138.1 d 


: 6.2 


Girardi et al. <20021 theoretical BCs (V filter) 


139.8 d 


: 5.3 


Girardi et al. <20021 theoretical BCs (K filter) 


138.8 d 


: 3.3 


Surface brightness-(6 - V) relation 


138 d 


: 19 


Surface brightness-r c ff relation (V filter) 


141.4 d 


: 9.9 


Surface brightness-r cff relation (K filter) 


139.1 d 


: 3.5 



errors. This is, however, a more relaxed constraint on the ef- 
fective temperature scale than that involved in finding distance 
using BCs. We note that our effective temperature uncertain- 
ties include contributions due to possible spectral peculiarity 
and systematic offset relative to the (inhomogeneous) effective 
temperatures used in the KTDS04 calibration. The uncertainty 
in distance could be reduced by further observations and es- 
timations of the effective temperatures of the two stars, using 
the same technique as for the stars used to calibrate the surface 
brightness relations. 

6. Conclusion 

The 'long' distance scale of the Pleiades is 132 + 3pc and is 
supported by main sequence fitting, the distance of the astro- 
metric binary Atlas, and by ground-based and Hubble Space 
Telescope parallax measurements. The 'short' distance is 120 + 
3 pc and is derived from parallaxes observed by the Hipparcos 
satellite. These results have been summarised in Table[S] It has 
been suggested that the two distances could be reconciled if the 
Pleiades cluster is metal-poor, but determinations of the atmo- 
spheric metal abundances of Pleiades F dwarfs suggest that the 
cluster has a solar iron abundance. 

We have studied the dEB HD 23642, a member of the 
Pleiades with an Hipparcos parallax determination, to calcu- 
late reliable absolute dimensions and uncertainties of the com- 
ponent stars. By comparing the radii of the components of 
HD 23642 to theoretical models we find that the metal and he- 
lium abundances are approximately solar, which removes the 
possibility that the 'long' and 'short' distance scales could be 
reconciled by adopting a low metal abundance or high helium 
abundance for the Pleiades. 

We have investigated the use of bolometric corrections 
(BCs) for determining the distances to eclipsing binaries, us- 



ing the empirical BC calibration of Code et al. ( 1976 ) and two 
sources of theoretically-calculated BCs. We find that the em- 
pirical and theoretical BCs give distances to HD 23642 in good 
agreement with each other. Distances determined using BCs 
for near-infrared filters are more precise and reliable due to a 
smaller dependence on interstellar reddening and metal abun- 
dance. 

We have presented a new, almost entirely empirical, tech- 
nique for determining the distance to dEBs composed of two 
components with effective temperatures between 10000K and 
and 3600 K, based on interferometrically-derived calibrations 
between effective temperature and surface brightness (Kervella 
et al. I2004> . This method does not explicitly require a light 
ratio for calculation of distance. Distances determined using 
the near-infrared JHKL calibrations are more precise as un- 
certainties in interstellar reddening and effective temperature 
are less important. Using this technique and Zf-filter photom- 
etry from 2MASS, we find that HD 23642 is at a distance of 
139.1 + 3.5 pc. This distance is consistent with the 'long' dis- 
tance scale of the Pleiades but in disagreement with the dis- 
tances to HD 23642 and to the Pleiades derived from Hipparcos 
parallax observations Table[S}- 

Further observations of HD 23642, to determine more ac- 
curate dimensions, would provide very precise metal abun- 
dance and effective temperature measurements for the compo- 
nent stars. This would reduce the uncertainty in its distance 
and allow further investigation of the system, which is itself an 
interesting object due to the metallic-lined nature of the sec- 
ondary star. Further infrared observations would also allow the 
use of entirely empirical surface brightness relations for the cal- 
culation of an accurate distance to HD 23642. 
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